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ABSTRACT 

Supernova remnants (SNRs) in the Galaxy are an important source of energy injection into 
the interstellar medium, and also of cosmic rays. Currently there are 294 known SNRs in the 
Galaxy, and their distribution with Galactocentric radius is of interest for various studies. Here 
I discuss some of the statistics of Galactic SNRs, including the observational selection effects 
that apply, and difficulties in obtaining distances for individual remnants from the ‘Z-D" 
relation. Comparison of the observed Galactic longitude distribution of a sample of bright 
Galactic SNRs - which are not strongly affected by selection effects - with those expected 
from models is used to constrain the Galactic distribution of SNRs. The best-fitting power- 
law/exponential model is more concentrated towards the Galactic centre than the widely used 
distribution obtained by Case & Bhattacharya (1998). 

Key words: supernova remnants - Galaxy: structure - ISM: structure 


1 INTRODUCTION 

There are currently 294 supernova remnants (SNRs) known in the 
Galaxy (Green 2014a). These are an important source of the injec¬ 
tion of energy and heavy elements into the interstellar medium, and 
are believed to be the site of acceleration of cosmic rays, at least up 
to 10*^ eV (e.g. Bell 2014). Consequently the distribution of SNRs 
in the Galaxy is of interest for a variety of studies of the Galaxy 
(e.g. Lee et al. 2011; Vladimirov et al. 2012; Kumar & Eichler 
2014; Galore etal. 2015). 

It is not straightforward to construct a Galactic distribution 
directly from properties of catalogued SNRs. This is not only be¬ 
cause distances are not available for most Galactic SNRs, but also 
because of the observational selection effects that apply to the cur¬ 
rent catalogue of SNRs. It is expected (Li et al. 2011) that most 
SNRs are from massive stars (i.e. types Il/Ib/Ic), which also pro¬ 
duce pulsars, and hence the distribution of SNRs will be closely 
related to that of star forming regions or pulsars. 

Some statistics of the current Galactic SNR catalogue, the se¬ 
lection effects that apply to their identification, and issues related to 
the derivation of distances to individual SNRs from the ‘Z-D’ rela¬ 
tion are discussed in Section 2. Section 3 presents constraints on the 
Galactic distribution of SNRs from comparison of the /-distribution 
of bright SNRs with various models. This includes comparison with 
the distribution from Case & Bhattacharya (1998), hereafter CB98 
(see also Case & Bhattacharya 1996), which has been widely used. 
The conclusions are summarised in Section 4. Preliminary results 
have been presented in Green (2012, 2014b), but here a more de¬ 
tailed analysis is made. 
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2 GALACTIC SNRS 
2.1 The catalogue of SNRs 

I have produced several catalogues of Galactic SNRs, with pub¬ 
lished versions in Green (1984, 1988, 1991, 1996a, 2004, 2009, 
2014a), and also various versions online since 1993'. Of the 
294 remnants in the most recent version of the catalogue, most 
have been first identified at radio wavelengths, or - if identified 
at other wavebands - have subsequently been detected at radio 
wavelengths. But there are 22 remnants in the catalogue that have 
not yet been detected at radio wavelengths, or not sufficiently well 
observed to provide an integrated radio flux density. At optical and 
X-ray wavelengths only about 40% and 30%, respectively, of the 
catalogued SNRs have been detected (which is not surprising, due 
to Galactic absorption that affects these wavelengths). Thus, it is 
selection effects at radio wavelengths that are dominant when con¬ 
sidering the completeness of the catalogue. Figure 1 shows a his¬ 
togram of the date of inclusion of a remnant in the catalogue (the 
larger number of entries in the first version of the catalogue, from 
1984, were largely taken from earlier Galactic SNR catalogues). 
This shows that major radio surveys have been the cause of signi¬ 
ficant increases in the identification of Galactic SNRs, notably: 

(i) the Effelsberg 2.7-GHz survey (Reich et al. 1990; Fiirst et al. 
1990), covering -2.6° < / < 240°, \b\ < 5°, with a resolution of 

4.3 arcmin; 

(ii) the MOST 843-GHz survey (Green et al. 1999), covering 
245° < / < 255°, \b\ < L5, with a resolution of x: 0.7 arcmin at 
best; 

' see: http://www.mrao.cam.ac.uk/surveys/snrs/ 
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Year of cataloging 

Figure 1. Histogram of the dates when Galactic SNRs were first catalogued. 



Angular size 6 / arcmin 

Figure 2. Histogram of the angular size, S, of the catalogued Galactic SNRs 
(there are also 10 large remnants with angular sizes greater than 2° which 
are not included in this plot). 


(iii) Brogan et al. (2006)’s survey of a small region of the 1st 
Galactic quadrant - 415 < I < 22°, \b\ < 1125 - at multiple radio 
wavelengths, with a resolution of « 0.6 arcmin at 327 MHz, which 
also used infrared observations to help discriminate between differ¬ 
ent types of sources. 

In the latest revision of the Galactic SNR catalogue, 21 new rem¬ 
nants were added to the catalogue. Of these 13 have an integrated 
radio flux density at 1 GHz in the catalogue, and hence a surface 
brightness, E can be calculated; all of these 13 remnants are faint, 
with El GHz < 8 X 10“^' W m“^ Hz“‘ sr“'. The other 8 newly cata¬ 
logued remnants do not yet have radio flux densities reported in the 
literature, as they were identified at optical or X-ray wavelengths 
(e.g. G38-7-1-3, G65-8-0-5, G66-0-0-0, G67-6-0-0, G67-6-)-0-9 
and G67-8-f0-5 were identified from optical observations by Sabin 
et al. 2013 - for these SNRs integrated radio flux densities have not 
been yet published). 


In order to derive directly the distribution of SNRs in the 
Galaxy both: (i) correction for the incompleteness of the current 
catalogue of SNRs due to observational selection effects, and (ii) 
knowledge of the distances to each known SNR is required. I dis¬ 
cuss below the selection effects that apply, at radio wavelengths, to 
the current SNR catalogue, and also some issues with the derivation 
of distances to individual SNRs using the ‘E-D’ relation. 


2.2 Selection effects 

As noted above, most Galactic SNRs are identified af radio 
wavelengths, and the selection effects that apply - as discussed 
previously, e.g. Green (1991, 2005) - are (i) intrinsically faint rem¬ 
nants (i.e. low surface brightness) are difficult to identify, and (ii) 
physically small but distant remnants are difficult to recognise as 
SNRs, due to their small angular sizes. 

I have previously argued (Green 2004) for a nominal surface 
brightness completeness limit of 10“^° W rcT^ Hz“* sr“*. (which 
is equivalent to 65 mjy per 1-arcmin circular beam). However, 
a higher limit probably applies close to the Galactic Centre (GC), 
where the background Galactic radio emission is brighter than else¬ 
where. The validity of this value for the approximate surface bright¬ 
ness completeness of the current SNR catalogue is discussed fur¬ 
ther below. 

In addition to the difficulty of identifying low-surface bright¬ 
ness SNRs, it is also difficult to identify small angular size rem¬ 
nants. For example, the Effelsberg 2.7-GHz survey has a resolution 
of 4.3 arcmin, and SNRs would have to be several times this angular 
size for their structure to be recognised. Reich et al. (1988) repor¬ 
ted 32 new supernova remnants in the first part of the area covered 
by the Effelsberg 2.7-GHz survey (-2';6 < / < 76°, \b\ < 5°, sub¬ 
sequently published in Reich et al. 1990). Of these 30 are resolved, 
with angular diameters of 16 arcmin or more, i.e. several times the 
resolution of the survey. The other two sources have small angular 
size, and were thought to be SNRs on the basis of other, higher res¬ 
olution targeted observations: (i) G54T-I-0-3, then thought to be a 
small, X! 1.5 arcmin ‘filled centre’ remnant (see Green 1985; Reich 
et al. 1985), but is now catalogued as a slightly larger 12 ar- 
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Table 1. SNRs in Effelsberg 2.7 GHz survey region with Zi ghz > 10 
W rrT^ Hz“' sr“' identified since 1992. 


name 

6 

/ arcmin 

S 1 GHz 

/Jy 

logiil GHz 

/ W nT^ Hz“* sr“* 

GO'3+O-O 

15x8 

22 

-19.56 

GlO-01 

8 

15 

-19.45 

G6-5-04 

18 

27 

-19.90 

G12-8-0-0 

3 

0.8 

-19.87 

G18-I-0-1 

8 

4.6 

-19.97 

G204+0-1 

8 

9.0 

-19.67 


cmin composite remnant, since a faint X-ray halo was detected by 
Bocchino et al. (2010), and (ii) G70-7-t-l-2, which was reported as 
a SNR by Reich et al. (1985), but this identification was not sup¬ 
ported by subsequent observations (e.g. Green 1986; Onello et al. 
1995; Cameron & Kulkarni 2007 and references therein). Given 
that SNRs are expected to have a continuous range of physical dia¬ 
meters, and are seen at a range of distances in the Galaxy, then a 
smooth distribution of angular sizes is expected. Figure 2 shows a 
histogram of the angular size, 6, of catalogued SNRs, which clearly 
shows a sharp decrease at small angular sizes (< 4 arcmin), due to 
the difficulty in identifying small SNRs. 

Although there are some small angular size remnants in the 
catalogue identified from high resolution observations - notably 
Gl-9-)-0-3 (e.g. Green et al. 2008; Reynolds et al. 2008; Borkowski 
et al. 2014), with an angular diameter of only 1.5 arcmin, many are 
missing. These missing young but distant SNRs will be on the far 
side of the Galaxy, and hence concentrated towards I = 0°,b = 0°, 
where confusion due to other Galactic sources along the line of 
sight makes their identification difficult. 

To justify the nominal surface brightness completeness limit 
of X! lO^^** W Flz^' sr“' is appropriate. Fig. 3 shows histograms 
of: (a) 272 of the 294 catalogued SNRs which have an integrated 
flux density, and hence a surface brightness at 1 GFIz; (b) 171 SNRs 
in the Elfelsberg 2.7 GFlz survey region; (c) 38 SNRs in the Effels¬ 
berg 2.7 GHz survey region entered into my SNR catalogue since 
1992, i.e. those not identified from the Effelsberg 2.7 GHz survey. 

Eigure 3(a) shows that the majority of SNRs identified are 
fainter than this nominal surface brightness limit, with 69 ‘bright’ 
remnants above the limit. The fainter remnants are more easily de¬ 
tected in regions where the Galactic background emission is fainter, 
i.e. away from b = 0° and 1 = 0°, which is illustrated in Eigs 4 and 
5. Eigure 4 shows the distribution in Galactic coordinates of (a) all 
remnants, and (b) the 69 ‘bright’ remnants. The distribution of all 
remnants, without taking the surface brightness effect into account, 
is very much broader in both coordinates than that of the bright 
remnants (the r.m.s. deviations from the Galactic Centre are 42° and 
2°4 in / and b respectively for all SNRs, with 36° and F3 for the 69 
‘bright’ SNRs). Eigure 5 shows the surface brightness and Galactic 
longitude of the 272 catalogued SNRs with radio flux densities. In 
the Galactic anti-centre (i.e. the 2nd and 3rd Galactic quadrants) - 
where the Galactic background is low - there is a higher proportion 
of low surface brightness SNRs. 

Eigure 3(b) and (c) show that the majority of SNRs identi¬ 
fied in the Effelsberg 2.7 GHz survey region from subsequent other 
observations are fainter than the nominal surface brightness limit. 
However, there are 6 SNRs above the surface brightness limit, 
which have subsequently been identified, see Table 1 (some of these 
had been suggested as possible SNRs earlier). These 6 remnants 
are all close to the Galactic centre, or are small angular size, so 


a) 



b) 



c) 



Figure 3. Histograms of surface brightness at 1 GHz for catalogued SNRs 
with a radio flux density at 1 GHz for: (a) SNRs; (b) SNRs in the region 
covered by the Effelsberg 2.7-GHz survey (i.e. -2° < 240°, \b\ < 5°), (c) 
SNRs in the Effelsberg suiwey region identified since 1992. 
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Figure 4. The distribution of SNRs in Galactic coordinates of all 294 catalogued SNRs for: (a) all SNRs, (b) the 69 ‘bright’ SNRs with a surface brightness 
above 10““^ W Hz“^ sr“^. Note that the latitude scale is exaggerated. 
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Figure 5. Radio surface brightness at 1 GHz against Galactic longitude for the 272 catalogued SNRs with radio flux densities. The nominal surface brightness 
limit of X\ GHz = 10“^® W m~^ Hz“^ sr“^ is shown. 
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a) 



Figure 6. Histogram of the distribution of Galactic SNRs (left scale) and 
the cumulative fraction (right scale) with Galactic longitude, for (a) all 294 
catalogued SNRs, and (b) the 69 ‘bright’ SNRs. 


would not have been well resolved in the EfFelsberg 2.7 GHz sur¬ 
vey, which is necessary to be recognised as a SNR. 

A surface brightness completeness limit of Zi ghz ~ 
10“^° W Hz ' sr“' is also consistent with observations made 
by Xu et al. (2013). Their observations covered 66“’ < / < 90°, 
\b\ < 4° at 5 GHz, and they used other radio and infra-red observa¬ 
tions to separate non-thermal and thermal components. This multi¬ 
wavelength approach should allow fainter remnants to be recog¬ 
nised than is possible using a single observation frequency. They 
concluded that there were no large, i.e. sufficiently resolved, rem¬ 
nants with El GHz ^ 0-37 X 10“^° W m“^ Hz“* sr“* in the region they 
observed. 

Figure 6 shows the distribution of SNRs with Galactic longit¬ 
ude. For all SNRs, Fig. 6(a) shows a clear asymmetry on different 
sides of ( = 0°, with more SNRs identified in the 1st and 2nd quad¬ 
rants. This is because (i) the 1st and 2nd quadrants are accessible to 
a large number of northern hemisphere telescopes, and (ii) they in¬ 
clude the deep, multiwavelength survey of 4°5 < / < 22°, \b\ < l‘;25 
by Brogan et al. (2006), as discussed in Section 2. The sample 
of 69 ‘bright’ SNRs, Fig. 6(b), shows a much smaller asymmetry 
close to / = 0°, but this is not statistically significant (4 in the re¬ 
gion 350° < / < 0°, and 8 - including GO-O-l-0-0 in the region 
0° < Z < 10°). By Galactic quadrants, the numbers of bright SNRs 
are: 35 in the 1st, 3 in the 2nd, 2 in the 3rd, and 29 in the 4th. 
These, with Poisson errors, do not show any asymmetry between 
either side of the Galactic Centre. This implies that the nominal 
completeness surface brightness cutoff, although derived from the 
EfFelsberg 2.7-GHz - which misses almost all of the 4th quadrant - 
is applicable there. Figure 6 again shows that the distribution of the 
‘bright’ SNRs is more localised towards Z = 0° than the distribution 
of all catalogued SNRs. 


2.3 The ‘E-D’ relation 

Although distance determinations are available for some Galactic 
SNRs - e.g. Pavlovic et al. 2013 give a compilation of distances 


to 60 SNRs from the literature - the E-D relation has been used 
for some time (e.g. Ilovaisky & Lequeux 1972; Clark & Caswell 
1976) to determine distances for other remnants. This is based on 
the fact that, for remnants with known distances, the observed sur¬ 
face brightness (27) is larger for SNRs with smaller physical dia¬ 
meters (D). This correlation is parameterised as 

27 = AD", (1) 

i.e. a straight line in the log D- log 27 plane, with A and n determ¬ 
ined from the properties of SNRs with known distances. The ob¬ 
served surface brightness, 27 - which is distance-independent - for 
a remnant without a distance determination can then be used to de¬ 
termine its physical diameter, D. Hence via its observed angular 
size, 9, its distance d = O/fZ can be determined. However - as pre¬ 
viously discussed (Green 1991, 2005) - there is a large scatter in 
the observed E-D distribution of SNRs, about an order of mag¬ 
nitude in D for a given 27. In addition, small and/or faint SNRs are 
more likely to have been missed in current surveys, and be missing 
from the current Galactic SNR catalogue, due to the observational 
selection effects discussed in 2.2. So, the true range of diameters 
for a given surface brightness may extend to lower diameters. (On 
the other hand, for a given surface brightness the upper limit to the 
range of diameters is not affected by selection effects, at least down 
to the nominal surface brightness limit of current catalogues.) Al¬ 
though the upper-right boundary of the observed E-D distribution 
of SNRs with known distances is useful to provide an upper limit on 
the diameter of an individual SNR, distances for individual SNRs 
derived from the E-D relation are imprecise. 

Leaving aside issues with selection effects, if a E-D relation is 
to be used statistically to derive distances to individual SNRs from 
their surface brightness, care has to be taken to use the appropriate 
form of regression (Green 2005). Since the E-D relation is used to 
derive D from 27 then a regression that minimises square deviations 
in logZ) should be used, not a regression than minimises square 
deviations in log 27. For distributions with a large spread - as is the 
case for SNRs with known distances - different regressions give 
significantly different fits (see Isobe et al. 1990; Feigelson & Babu 
1992 and Feigelson & Babu 2011 for discussion of different forms 
of least square regressions). 

CB98 derived a E-D relation with 

27 oc O-2 f>4±0.26^ 

using 37 Galactic ‘shell’ SNRs with known distances (or 27 oc 
£,-2.38+0.26 if (=Glll-7-2-l) is excluded). CB98 comment 

that these E-D slopes are significantly flatter than those derived 
in previous studies, e.g. Milne (1979) obtained 27 oc D^^^. It is 
evident, however, that CB98 minimised square deviations in log 27, 
whereas Milne (1979) minimised square deviations in logD, so 
these E-D slopes are not directly comparable. Re-fitting the SNRs 
with known distances used by CB98, minimising deviations in 
log£) rather than log 27 produces much steeper 27- Z9 slopes, which 
are in good agreement with earlier results (e.g. those of Milne). For 
all 37 ‘shell’ SNRs 


27oc (3) 

(or 27 oc /)-3-37 ±o. 35 jf ^ excluded). This difference in the 
slope of the derived E-D relation depending on the form of re¬ 
gression used has an important consequence for the Galactic SNR 
distribution derived by CB98. There is a systematic bias in the dia¬ 
meters and hence distances derived for SNRs, with fainter/brighter 
SNRs having distances that are too large/small respectively. Since 
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Figure 7. Histogram of the distribution of Galactic SNRs (left scale) and the 
cumulative fraction - solid black line - (right scale) with Galactic longitude 
for the 69 ‘bright’ SNRs. The dashed blue line is the cumulative fraction for 
the best-fitting power-low/exponential model (with a = 1.09 and /? = 3.87). 


Figure 8. As Fig. 7, but the dashed blue line is the cumulative fraction for 
the best-fitting power-low/exponential model from CB98 (with a = 2.00 
andy3 = 3.53). 


there are more faint SNRs than bright SNRs, then the derived distri¬ 
bution of SNRs in the Galaxy will be systematically spread out too. 
In E-D studies in the literature it is not always clear what form of 
regression was used. Re-analysing the published lists of SNRs with 
known distances implies that others (e.g. Gobel et al. 1981; Huang 
& Thaddeus 1985; Arbutina & Urosevic 2005; Stupar et al. 2007) 
have-, like CB98, used regressions which minimise deviations on 
log 27 rather than log D. 

Pavlovic et al. (2013) discussed the fitting of a E-D relation to 
Galactic SNRs with known distances, and concluded from Monte 
Carlo simulations that ‘orthogonal regression’ gives the best res¬ 
ult. This conclusion is to be expected, given that in their Monte 
Carlo simulations Pavlovic et al. assumes errors in logD and log 27 
of same magnitude, which is not realistic. In practice, the errors 
in surface brightness - which depends on the uncertainties in flux 
density and angular size - are likely to be smaller than the errors 
in distances (as Pavlovic et al. say in their introduction). When the 
errors in logZl and log 27 are different, then the orthogonal fitting 
does not produce the correct result. However, if the E-D relation is 
to be used to derive distances to individual SNRs, then a regression 
minimising deviations in logD should be used to obtain the best 
result. 


3 THE GALACTIC RADIAL DISTRIBUTION OL SNRS 

Given the limitations of the E-D relation and the selection effects 
that apply to the identification of SNRs, rather than deriving the 
distribution of SNRs in the Galaxy directly - as done by CB98 - 
an alternative approach is to consider the /-distribution of sample 
SNRs. This can be compared with the expected /-distribution from 
various models. 

Here I choose the sample of 69 ‘bright’ SNRs, above the nom¬ 
inal surface brightness limit of 27| ghz = 10“^° W m“^ Hz“’ sr“*, 
where the current catalogue of Galactic SNRs is thought to be 
nearly complete. This method has the advantage of avoiding all of 
the uncertainties from the E-D relation, but has the disadvantage 
of only using a fraction of the total number of catalogued Galactic 
SNRs. I include checks on how dependent the results are on a 
choice of surface brightness limit, or exclusion of regions close to 
/ = 0°. As noted above, a single surface brightness limit is not ap¬ 
propriate at all Galactic latitudes, and some SNRs above this nom¬ 
inal limit may well be missed in the brightest part of the Galactic 


^ This is contrary to the statement in Pavlovic et al. (2013) that all previous 
E-D studies used regression minimising deviations in log 27. 



Figure 9. Normalised surface density of SNRs with Galactocentric radius 
for power-low/exponential models: (i) solid line, the best fit to /-distribution 
of 69 bright SNRs, and (ii) dashed line, from CB98. 


plane, i.e. at latitudes near 1 = 0°. Moreover, the difficulty of identi¬ 
fying small angular size remnants is also more of an issue nearer 
1 = 0°. Nevertheless, this sample is not as strongly affected by se¬ 
lection effects as the complete catalogue. 

Here I use a power law/exponential model for the Galactic 
surface density of SNRs with Galactocentric radius, R, with 


oc 




{R-Re) \ 

Re )’ 


(4) 


where Rq = 8.5 kpc (i.e. a cylindrically symmetrical distribution 
about the Galactic Centre). Given the limited number of SNRs in 
the sample of ‘bright’ remnants, no attempt is made here to con¬ 
strain the form of the distribution perpendicular to the Galactic 
plane. This model for the surface density of SNRs tends to a zero 
density towards the Galactic Centre, which better matches the dis¬ 
tributions derived for pulsars and star formation in the Galaxy (e.g. 
Johnston 1994; Bronfman et al. 2000; Paladini et al. 2004; Yusifov 
& Kiifuk 2004) than, for example, a simple Gaussian distribution 
(which I used in Green 1996a). This power law/exponential is one 
of the two models used by CB98, who obtained best-fitting para¬ 
meters of ff = 2.00 ± 0.67,= 3.53 ± 0.77. 

The best-fitting power-law/exponential model to the 69 bright 
remnants, is shown in Fig. 7, which has a = 1.09 and p = 3.87. 
The statistic used for the fitting was the least sum of squares of the 
differences between the observed and model cumulative fractions, 
summed over all latitude bins, i.e. if the observed and model cu¬ 
mulative fractions are /o(/) and /m(/) respectively, for the /th bin, 
the minimum of 2,(/o - fmf- (The fitting was made with a model 
distribution of one million SNRs, out to a Galactocentric radius 
of 50 kpc, using bins of 10“’ in /.) The sum of squares misfit for 
the a = 1.09 and p = 3.87 bets fitting model is 0.0127. Figure 8 
shows the power-law/exponential model derived by CB98, which 
is a poor fit to the observed distribution of the 69 bright SNRs, with 
a sum of squares misfit of 0.0841. This clearly represents a broader 
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Table 2. Best-fitting power-law/exponential models for different samples of Galactic SNRs, and other power-law/exponential models for comparison. 



Zl GHz cut-off 

/-range 

number 

fit parameters 

proportion inside 


/Wm-2 Hz-' sr-' 


of SNRs 

a 

P 

Solar Circle 

bright sample 

C 

1 

O 

all 

69 

1.09 

3.87 

73% 

brighter sample 

2 X 10-2° 

all 

44 

0.51 

2.91 

70% 

fainter sample 

5 X 10-2' 

all 

103 

1.49 

4.60 

77% 

omit near GC sample 

O 

1 

O 

|/| > 10° 

57 

0.0 

2.76 

77% 

CB98 best fit 

o 

1 

O 

all 

69 

2.00 

3.53 

49% 

CB98 a fixed 

o 

1 

O 

all 

69 

2.00 

5.11 

76% 

CB98 j8 fixed 

o 

1 

O 

all 

69 

0.85 

3.53 

73% 


Notes: (a) for the first four rows, the a and /? values are for the best ht; (b) the hfth row uses the a and /? 
best-fitting values from CB98; (c) the final two entries use one or other of the bestOfitting values for a and 
j8 from CB98 and varies the other for a best fit. 


distribution (e.g. a flatter slope of the cumulative distribution near 
1 = 0° than what is observed for the ‘bright’ remnants) than that 
shown in Fig. 7. This is not surprising given the systematic bias in 
the distances derived from the Z-D relation used by CB98 which 
was noted in Section 2.3. The parameters for the best-fitting power- 
law/exponential model for the sample of 69 bright SNRs, and that 
from CB98 are shown in Table 2. The broader CB98 distribution is 
also evident in Fig. 9, which shows the surface density distribution 
of SNRs with Galactocentric radius for CB98’s model, and the best 
fit to the 69 bright remnants. For CB98’s best-fitting model, only 
49 per cent of SNRs are inside the Solar Circle, compared with 
73 per cent for the best fit to the 69 bright remnants. 

The parameters of the best-fitting model are not well defined, 
as there is a strong degeneracy between the parameters, e.g. it is 
possible to obtain almost as good fits with, a = 2.00, J3 = 5.11 or 
a = 0.85,= 3.53 (i.e. fixing one parameter to the value obtained 
by CB98, and varying the other), which have sum of squares misfits 
of 0.0135 and 0.0128 respectively. The distribution with a = 2.00, 
j8 = 5.10 peaks at a slightly larger radius from GC than the best¬ 
fitting distribution with a = 1.10 and jj = 3.90, but has a similar 
fraction of SNRs, 76 per cent, inside the Solar Circle. Conversely 
the distribution with a = 2.00, /3 = 5.10 peaks closer to the GC, but 
again has a similar fraction of SNRs, 73 per cent, inside the Solar 
Circle. 

Also, the best-fitting model depends on (i) the accuracy of the 
value of the surface brightness cut-off chosen to define the sample 
of ‘bright’ SNRs, and also (ii) the assumption that the remaining 
selection effects, which are more important close to f = 0° are not 
important. To investigate these I have done the analysis with differ¬ 
ent Zi GHz cut-off values (lower and higher by a factor of two), and 
excluding SNRs within 10° of / = 0'’. Varying the surface bright¬ 
ness by a factor of two - the ‘brighter’ and ‘fainter’ samples in 
Table 2 - does give rather different parameters, and hence different 
radial distribution, but the fraction of SNRs within the Solar Circle 
does not change strongly. When the region |/| < 10° is excluded 
- the ‘omit near GC’ sample in Table 2 - the best-fitting model 
is a pure exponential (but note that this fitting does not depend on 
R < 1.5 kpc). This suggests there are indeed residual selection ef¬ 
fects near 1 = 0°. Nevertheless, the main result, that the distribu¬ 
tion obtained by CB98 is too broad compared with bright SNRs 
still holds. The conclusion is strengthened by the fact that any re¬ 
maining incompleteness of the sample of 69 bright SNRs will be 
concentrate towards 1 = 0°. 

Finally, it should be noted that the constraints presented here 
apply to the distribution of observed SNRs. If the observability of 
SNRs above the radio surface brightness limit depends on a prop¬ 


erty - e.g. ambient density - which varies with Galactocentric ra¬ 
dius, then this will mean the distribution of SNRs is not the same 
as their parent SNe. 


4 CONCLUSIONS 

Here I have discussed some of the properties of the most recent 
catalogue of Galactic SNRs, which contains 294, particularly (i) 
the selection effects that apply to the completeness of the catalogue, 
and (ii) issues with using the ‘27-/9’ relation to derive distances to 
individual remnants. By comparison of the distribution in Galactic 
longitude of a sample of 69 ‘bright’ SNRs - which are not strongly 
affected by selection effects - with that expected from models, con¬ 
straints are placed on the distribution of the SNRs with galacto¬ 
centric radius (using a power-law/exponential model). It is shown 
that the widely used distribution derived by Case & Bhattacharya 
(1998) is too broad. 
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